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Abstract. We discuss the first results of our ongoing pro- 
gram of Target of Opportunity observations of gamma-ray 
bursts (GRBs) using the SCUBA instrument on the James 
Clerk Maxwell Telescope. We present the results for GRB 
970508, 971214, 980326, 980329, 980519, and 980703. 

Our most important result to date is the detection 
of a fading counterpart to GRB 980329 at 850 fxm. Al- 
though it proved to be difficult to find the infrared coun- 
terpart to this burst, the sub- millimeter flux was relatively 
bright. This indicates that intrinsically the brightness of 
this counterpart was very similar to GRB 970508. The 
radio through sub-millimeter spectrum of GRB 980329 is 
well fit by a power law with index a = +0.9. However, 
we cannot exclude a ^ 1/3 power law attenuated by syn- 
chrotron self-absorption. An a ~ +1 VLA-SCUBA power 
law spectrum is definitely ruled out for GRB 980703, and 
possibly also for GRB 980519. 

We cannot rule out that part of the sub- millimeter flux 
from GRB 980329 comes from a dusty star-forming galaxy 
at high redshift, such as the ones recently discovered by 
SCUBA. Any quiescent dust contribution will be much 
larger at sub- millimeter than at radio wavelengths. Both a 
high redshift and large dust extinction would help explain 
the reddening of the counterpart to GRB 980329, and a 
redshift of z ~ 5 has been suggested. The large intensity of 
this burst might then indicate that beaming is important. 

Key words: gamma rays: bursts - radio continuum: gen- 
eral - infrared: general 



1. Introduction 

An important reason for the long-standing mysteries sur- 
rounding the gamma-ray burst (GRB) sources has been 
the lack of prompt accurate locations to look for quies- 
cent or fading counterparts. This situation changed dra- 
matically in 1997 with the rapid and accurate location 
of X-ray emission detected during GRBs using the wide- 
field camera on the Satellite per Astronomia X (Bep- 
poSAX), the All Sky Monitor on the Rossi X-ray Timing 
Explorer (RXTE), and from the detection of X-ray after- 
glow in scans of BATSE GRB error boxes with the PCA 
on RXTE. 

These localized X-ray counterparts have led to intense 
multiwavelength campaigns. Optical transients have been 
found to some but not all of the bursts with fading X-ray 
counterparts. The reason optical transients have not been 
discovered for the other bursts with X-ray counterparts 
(some of which are much brighter during the gamma-ray 
burst itself) is that the searches have not always been sen- 
sitive enough, particularly given the problem of absorption 
local to the source (Groot et al. 



1998a) 



In spite of detailed searches, radio emission has only 
been detected for a few bursts. The first detection of a 



variable radio sourc e was GRB 970 508 (F rail et al. p. 997 ; 
Taylor et al. [19971 Galama et al. |l998a|). This ~ mJy 



source had a rising spectral index, and its rapid radio 
variability was most likely caused by interstellar scintilla- 
tion. The combined observations of GRB 970508 showed 
that the peak of the spectrum was in the sub-millimeter 



region (Galama et al. 1998b) with the radio, millimeter, 



and optical emission peaking days to weeks after the burst 
(Frail et al. |1997|; G alama et al |1998a| ; Gruendl et al 
Bremer et al. 11998; Pian et al. 11998; Pcdersen et al. 
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Castro-Tirado et al. 1998). Even before any counterparts 
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were found, two completely separate classes of models had 
suggested this would be the case: (1) cosmological fire- 
ball models (e.g. Paczyhski & Rhoads 1993; Katz 1994; 
Meszaros & Rees 1997), and (2) Compton scattering mod- 
els (e.g. Liang et al. 1997; Liang 1997). 

For GRB 970508 there are at least two breaks in the 



spectrum between 10 11 and 10 14 Hz (Galama et al. 1998b) 
making this a crucial region for the models. To obtain a 
complete picture of the nature of the burst counterparts, 
it is clear that one needs to cover the entire spectrum, and 
sub-millimctcr observations with a ~ mJy sensitivity are 
needed. This is particularly important since the optical 
emission can be suppressed by local absorption, and the 
radio emission can be self absorbed as well as scrambled 



by interstellar scintillation (Walker 1998 ), while "clean" 
observations ought to be possible at sub-millimeter wave- 
lengths. 

In this paper we discuss our ongoing program of Target 
of Opportunity observations using SCUBA on the James 
Clerk Maxwell Telescope. In §2 we discuss some of the 
technical features of SCUBA that make it well suited for 
performing counterpart searches. In §3 we present the re- 
sults of our observations to date on GRBs 970508 (which 
was a limited trial run), 971214, 980326, 980329, 980519, 
and 980703. In §4 we give a brief discussion. 

2. SCUBA Details 

SCUBA is the new sub-millimeter continuum instrument 
for the James Clerk Maxwell Telescope on Mauna Kea, 



Hawaii (for a review see Holland et al. 1998). It uses two 
arrays of bolometers to simultaneously observe the same 
region of sky, ~ 2.3' in diameter. The arrays are optimized 
for operations at 850 and 450 /im. Fully sampled maps of 
the 2.3' region can be made by "jiggling" the array. Thus 
this mode is appropriate for mapping the better localized 
GRB error boxes as well as those of the X-ray transients. 
This mode can also be used to look for extended quiescent 
counterparts. 

Deeper photometry can be performed using just the 
central pixel of these arrays. There are also dedicated pho- 
tometry pixels for 1100, 1350, and 2000 /im observations 
(that cannot be used at the same time as the arrays). This 
photometry mode is appropriate for well localized radio or 
optical transients. 

Scan mapping of larger GRB error boxes using the 
450:850 filters is also possible. While this mode is ideal for 
mapping the long thin GRB error boxes that are obtained 
using triangulation between satellites, e.g. 5' x 0.5' (Hurley 
et al. 1997 ), the sensitivity is greatly reduced. 

Given the rapid dissemination of candidate optical and 
radio transients, the photometry mode is the one that we 
use most often. In principle, the most sensitive measure- 
ments can be made at 1350 /Ltm, though the 850 fj,m array 
has an advantage because the multiple bolometers permit 
a good sky noise subtraction. In photometry mode, for an 



integration time of 2 hours, we would expect to achieve 
an rms ~ 1 mJy at 1350 /im, ~ 1.5 mJy at 850 /xm, 
and ^5 — 20 mJy at 450 /im. The sensitivities depend 
significantly on the weather, particularly at the shorter 
wavelengths. The jiggle maps give an rms that is a factor 
^2 — 3 times higher. 



3. Results of SCUBA Observations 

3.1. GRB 970508 

As a limited trial run, a 30 minute SCUBA observation 
of GRB 970508 was made on 1997 May 26 using the 1350 
fj,m photometry pixel. The weather conditions were very 
poor. No source was detected with an rms ~ 10 mJy. This 
result is consistent with the other millimeter observations 
of GRB 970508 (Gruendl et al. 
Shepherd et al. |1998|). 



199S: Bremer et al. 



1998 



3.2. GRB 971214 

Our preliminary SCUBA results on G RB 971214 were 
originally reported in Smith et al. (1997 ) . 

The BeppoSAX GRB Monito r was triggered on 1997 
December 14.97 UT (Hcisc et al. |1997[ ). A previously un- 
known fading X-ray source (1SAX J1156. 4+6513) was 
found inside the burst error circle (Antonclli et al. 1997). 
Consistent with this X-ray source an optical transient 
was found (e.g. Halpern et al 



1998; Kulkarni et al. 1998 



Gorosabel et al. 1998: Diercks 



Ramaprakash et al. 1998 
et al. 1998). A possible quiescent host to the transient was 
found with a redshift z = 3.42 (Kulkarni et al. |1998| ). No 
ra dio co unterpart has so far been seen (Ramaprakash et 
al. |1998D . 

We began our series of SCUBA observations on UT 
1997 December 16, before the optical transient was re- 
ported, and when the error box of 1SAX Jl 156.4+6513 
had a radius ~ 1' (L. Piro, private communication). We 
performed a 450:850 jiggle map of the whole error box. 
Fig. [j] shows the whole 850 //m map, which illustrates a 
typical SCUBA jiggle map. At 850 /im, the rms was 3 mJy 
in the central region of the map and the beam size was 
14.7". The optical transient was near the edge of this map 
at RA(J2000) = 11:56:26.4, DEC(J2000) = +65:12:00.5, 
where the rms is approximately 5 mJy. No sources were 
detected by SCUBA anywhere in the map. 

The remainder of our SCUBA observations were per- 
formed on the optical transient using the photometric 
mode with the 450:850 arrays. The results are given in 
Table [l]. We did not detect a sub-millimeter continuum 
source at the location of the optical transient. Combining 
all our photometric observations gives an rms of 1.0 mJy 
at 850 /im. 

The decay of the optical flux followed a power law 
oc t~ s with slope 5 ~ 1.2. In the simple adiabatic pis- 
ton model, a fireball produced by a one time impulsive 
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Fig. 1. SCUBA 850 /im jiggle map of the ~ V radius error 
box to 1SAX J1156.4+6513 on UT 1997 December 16. The 
beam size is 14.7". 



injection of energy in which only the forward blast wave 
efficiently accelerates particles predicts a power law spec- 
with energy spectral index (5 — 25/3 



trum S v oc , 

(Wijers, Rees, & Meszaros |1997). For GRB 971214, this 
would imply [3 — 0.8. This is in general agreement with 
the optical to X-ray slope. Extrapolating this power law 
gives a flux density 5s 1 mJy at 850/im on December 17. 
However, the optical spectrum alone has a much steeper 
spectrum, indicating that there is significant extinction 



loc al to the source (Halpern et al. 1998; Ramaprakash et 
al. 1998| ). The required correction is several magnitudes in 
the I-band, which similarly raises the prediction at 850/im. 
Our SCUBA limits then imply that there is a break be- 
tween the optical and sub-millimeter bands, as was found 
in GRB 970508. Infrared observations of GRB 971214 sug- 
gest that this break was at ~ 1 /im in the first few hours 
after the burst (Ramaprakash et al. 1998; Gorosabel et al. 



199? 



3.3. GRB 980326 

The BeppoSAX GRB Monitor was triggered on 1998 



March 26.888 UT (Celidonio et al. |1998| ). Although Bep- 
poSAX was unable to make an observation with the Nar- 
row Field Instruments and RXTE did not see any X-ray 
emission from the GRB error box (Marshall & Takeshima 
1998j ), a can didate optical transient was found (Groot 



et al. 1998b). This transient was notable for its unusu- 



Table 1. SCUBA observations of the optical transient to 
GRB 971214. There is a ~ 10% systematic uncertainty in 
all the flux densities. 



Observing Time 


850 ^m rms 


(UT 1997) 


(mJy) 


Dec 16 


5.0 


Dec 17 


1.4 


Dec 19 


1.9 


Dec 22 


1.3 



2.10 ± 0.13 (Eichelberger et al. [199$ Groot et al. 

25.5±0.5 



1998b). A constant underlying source with R, 
was found (Grossan et al. 1998| ; Djorgovski et al. 1998a; 



Groot et al. |1998bp . This burst was also interesting in that 
the gamma-ray spectrum during the burst was quite soft. 

We used SCUBA to make a short photometry observa- 
tion of the optical transient to GRB 980326 on 1998 March 
29. The source was not detected, with rms 2.5 mJy at 850 
/im and 25 mJy at 450 /zm. Since there was no report of 
a radio counterpart, and the much more interesting GRB 
980329 occurred at this time, we did not try to make any 
further observations of GRB 980326. 

On 1998 March 29 the counterpart had R c = 24.5, 
and the optical spectrum was poorly determined, with 
/3 = 0.66 ± 0.70 (Groot et al. |l998b| ). Extrapolating with 
(3 = 0.66 would give a flux density of 0.06 mJy at 850 /im, 
while using the la value of /3 = 1.36 gives 8.8 mJy at 850 
/im. Possible distortions to the optical spectrum from the 
underlying quiescent source, and the difficulty in deter- 
mining the extinction corrections add to the uncertainty 
in the counterpart spectrum. Thus we are currently un- 
able to make any statements about breaks in the optical 
to sub- millimeter spectrum for GRB 980326. 

3.4. GRB 980329 

Our preliminary SCUBA results on GRB 980329 were 



originally reported in Smith & Tilanus (1997) 



The BeppoSAX GRB Monitor was triggered on 1998 



March 29.156 UT (Frontera et al. |1998| ). This was the 
brightest burst that had been seen simultaneously by the 
BeppoSAX Wide Field Camera, with a peak flux ~ 6 
Crab in the 2-26 keV band. A fading X-ray source 1SAX 
J0702. 6+3850 was found using the BeppoSAX Narrow 



Field Instruments (in't Zand et al. 1998). Inside this X-ray 
error box, a variable radio source VLA J070238. 0+385044 
was found that was similar to GRB 970508 (Taylor et al. 
|l998a " " 



1998b) 



It was not until after the variable radio 
source was discovered that infrared observations found a 
fading c ounter part (Klose et al. 1998| ; Palazzi et al. 1998 



ally rapid optical fade, with a power law decay index 



Metzger 1998 ): this indicated that the optic al ex tinction 
was significant for this source (Larkin et al. 199S; Taylor 
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Fig. 2. Evolution of the 850 /tm flux seen by SCUBA for 
the counterpart VLA J070238.0+385044 to GRB 980329. 
The two plots show the same data plotted on linear-linear 
and log-log axes. The curves plot flux density oc t~ m where 
t is the time since the burst and m = 1 (solid), m = 2 
(short dashed), and m — 3 (long dashed). 



et al. 1998b). A possible host ga laxy was found at this 
location (Djorgovski et al. 1998b ). 

Starting on 1998 April 5, we made a series of observa- 
tions of VLA J070238.0+385044 using SCUBA. On April 
5.2 UT, we detected the source at 850 /mi with a flux den- 
sity of 5 ± 1.5 mjy. This source was confirmed on April 
6.2 with a flux density of 4± 1.2 mjy, resulting in an aver- 
age of 4.5 ± 1 mjy over the two days. The source was not 
detected at 450 /tin, with an rms of 10.0 mjy averaged 
over these two days. The 850 /im source was present in 
all our separate integrations, making us confident that it 
was real. A hint of a fading trend was confirmed by ob- 
servations on April 7.2, when the 850 /mi flux density was 
2 ± 0.8 mjy. Observations on April 8 gave 1.5 ± 1.0 mjy 
at 850 /xm, with no detection at 1350 /tm (the rms was 1.2 
mjy). Finally, the signal was 1.4 ± 1.0 mjy at 850 /tm on 
April 11. 

Assuming the sub-millimeter fluxes are due to the 
burst counterpart, they should represent "clean" measures 
of its intensity, unaffected by scintillation and extinction. 
Although the optical emission was significantly reduced 
in this burst, the radio and sub- millimeter observations 
show that intrinsically the brightness of this counterpart 
was s imilar to GRB 970508 (e.g. see Figure 2 of Palazzi 
et al. |1998| ). 



Fig. 3. The radio through sub-millimeter spectrum of 
GRB 980329. The time-averaged VLA an d OVR O points 
are taken from Figure 2 of Taylor et al. ( |l998b| ). The 3cr 
SCUBA upper limit at 450 /tm has been included. The 
solid curve is a power law oc v a with a = +0.9. The 
dashed curve is a u 1 / 3 power law attenuated by a syn- 
chrotron self- absorption with Vq = 13 GHz at r = 1. 
The dotted curve adds a v 3 power law to the synchrotron 
curve. 



Fig. I plots the evolution of the 850 /tm SCUBA flux. 
For a power law decay with the flux density oc t~ m where 
t is the time since the burst, the best fit power law index is 
m = 2.8. However, m is not tightly constrained: the 90% 
confidence interval ism = 0.9tor7i = 5.6. 

Fig. | adds the SCUBA results to the VLA-OVRO 



results presented in Figure 2 of Taylor et al. (1998b). Be- 
cause of the averaging of the rapidly varying radio data 
over several days, some caution is required in using this 
figure. Taylor et al. found that a power law S v oc v a with 
a = +0.9 gave the best fit to the VLA-OVRO data alone. 
The solid curve shows that this extends very well to the 
SCUBA results. The dashed curve shows that the popular 
v 1 / 3 power law (e.g. Katz 1994; Waxman 1997) attenu- 



ated by a synchrotron self-absorption component gives a 
much worse description of the longer wavelength emission 
for GRB 980329. However, the reduced xl= 2.6 for this 
fit, and the probability that a random set of data points 
would give a value of xt as large or larger than this is 
Q = 0.034. It is therefore not possible to exclude this 
model, and it will be important to study more bursts to 
determine whether there is a range of power law indices 
for a. 
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One way to slightly reduce the sub-millimeter flux of 
the counterpart in Fig. ^ would be if part of the flux comes 
from an underlying quiescent sub-millimeter source. There 
is a suggestion in Fig. ^ that the 850 /im flux is leveling 
off, but an instrument more sensitive than SCUBA will 
be required to see if such a quiescent source is present 
for GRB 980329. SCUBA has recently discovered several 
dusty star-forming galaxies at high redshifts (Hughes et 
al. |1998j; Barger et al. |1998|; Smail et al. |1998[). Both a 



flux of 1.6 mJy at 850 /jm at the time of our SCUBA 
observation. Unfortunately, our observation was made too 
late to determine if there was a break between the optical 
and millimeter bands in GRB 980519. 



3.6. GRB 980703 



BATSE trigger 6891 (Kippen et al. pL998| ) was also detected 
by the RXTE ASM on 1998 July 3.182 UT (Levine et al. 



high redshift and large dust extinction would help explain 
the reddening of the counterpart to GRB 980 329, and a 
redshift of z ~ 5 has been suggested (Fruchter 1998 ). The 
large intensity of this burst might then indicate that beam- 
ing is important. 



1998). BeppoSAX NFI observations of the RXTE ASM er- 



In dust models, one expects S v oc or u 4 (e.g. Dwck Bloom et al. |l998b| ; Djorgovski et al 



ror box located a fading X-ray source 1SAX J2359. 1+0835 
(Galama et al. 1998c, 1998oj ). A variable radio, infrared, 
and optical counterpart was found, as well as an under- 
lying galaxy with R = 22.6 and a redshift of 0.966 (e.g. 

1998cj ). 



& Werner 1981). Thus any quiescent dust contribution 



is very much larger at sub-millimeter than at radio wave- 
lengths. For illustrative purposes the dotted curve in Fig. || 
adds a quiescent S v oc f 3 component to the synchrotron 
curve. In this example, the quiescent flux density at 8.3 
GHz is only 0.02 /iJy. 

3.5. GRB 980519 

The BeppoSAX GRB Mo nitor was triggered on 1998 May 
19.514 UT (Muller et al. |1998|). A fading X-ray counter- 



part 1SAX J2322. 3+7716 was found, alt hough the X-ray 
decay was not monotonic (Nicastro et al. 1998 ) . A fading 
optical counterpart was also found, w hose p ower law decay 
was steep 8 ~ 2 (e.g. Jaunsen et al. 1998, Djorgovski et 



al. 1998c ). A very faint quiesce nt op tical source w as even - 
tually detected (Sokolov et al. |l998| , Bloom et al. |l998a| ). 
A variable radio source was found at the same location 



(Frail et al. |l998a[ ). 

This source was not in an ideal location for SCUBA 
observations, with the elevation never rising above 35°. 
Also, the weather conditions were very poor at this time, 
and the JCMT was locked into using a different instru- 
ment. This meant we were only able to make one SCUBA 
observation of GRB 980519 on 1998 May 27 The source 
was not detected, with flux density 0.9 ±1.8 mJy at 850 
/im and an rms of 80 mJy at 450 /im. 

The radio flux uncorrected for scintillation at the time 
of our SCUBA observation is not currently available. On 
May 22.3, the 8.3 GHz flux measured by the VLA was 0.1 
mJy. Extrapolating from this using a = +1/3 predicts a 
flux density of 0.35 mJy at 850 /im. On the other hand, 
extrapolating using a = +1 predicts a flux density of 4.3 
mJy at 850 /im. When the final radio results are available, 
it may be possible to determine whether this steeper slope 
is unacceptable for GRB 980519. 

It is believed that the optical extinction is small for 



this burst (Gal et al. 1998). Assuming the optical flux 
continued to decay with a power law of index 8 — 1.98, 
and extrapolating the optical spect rum a ssuming a power 
law index of (3 = 1.26 (Gal et al. 1998 ) would predict a 



SCUBA observed the radio counterpart on 1998 July 
10.5 UT. The source was not detected, with an rms of 
2.7 mJy at 1350 /im. A second observation was performed 
on 1998 July 15.6 UT. Again the source was not detected 
with an rms of 2.1 mJy at 850 /im and 24 mJy at 450 /im. 
Another observation was tried on July 16, but the weather 
conditions were too poor to produce any useful results. 

While the 4.86 GHz flux suffered from large variations, 
the 8.46 G Hz flux was steadier, with a mean of 0.94 mJy 
(Frail et al. 1998b ). Extrapolating from this using a = +1 
predicts flux densities of 39 and 25 mJy at 850 and 1350 
/im respectively. Our SCUBA results can definitely rule 
out this simple power law for GRB 980703. Extrapolating 
from the radio using a = +1/3 predicts flux densities of 
3.3 and 2.8 mJy at 850 and 1350 /im respectively. While we 
would expect to have detected 1— 2cr signals in our SCUBA 
observations, the lack of detections are not inconsistent 
with this model for GRB 980703. 

4. Discussion 

The sub-millimeter is an important band for GRB studies 
because it is where the emission peaks in some bursts in 
the days to weeks following the burst. The sub-millimeter 
emission is not affected by extinction local to the source 
or interstellar scintillation. We have shown that sub- 
millimeter observations are important to: 

— Determine the radio to sub-millimeter spectrum, and 
see if it agrees with the synchrotron models. 

— Discover whether there are breaks in the optical to 
sub- millimeter spectrum. 

— Determine the evolution of the sub-millimctcr flux. 

— Look for underlying quiescent sources that may be 
dusty star-forming galaxies at high redshifts. 

To obtain a detailed understanding of the GRB 
counterpart behaviors will require observations of many 
bursters. To this end our program of Target of Opportu- 
nity observations using SCUBA is ongoing. 
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